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A fully three-dimensional (3D) model of an electron beam 
focused by a periodic permanent magnet (PPM) stack has been 
developed. First, the simulation code MAFIA was used to 
model a PPM stack using the magnetostatic solver. The exact 
geometry of the magnetic focusing structure was modeled; 
thus, no approximations were made regarding the off-axis 
fields. The fields from the static solver were loaded into the 3D 
particle-in-cell (PIC) solver of MAFIA where fully 3D 
behavior of the beam was simulated in the magnetic focusing 
field. The PIC solver computes the time-integration of 
electromagnetic fields simultaneously with the time integration 
of the equations of motion of charged particles that move under 
the influence of those fields. Fields caused by those moving 
charges are also taken into account; thus, effects like space 
charge and magnetic forces between particles are fully 
simulated [']. The electron beam is simulated by a number of 
macro-particles. These macro-particles represent a given 
charge Q amounting to that of several million electrons in order 
to conserve computational time and memory. Particle motion 
is unrestricted, so particle trajectories can cross paths and move 
in three dimensions under the influence of 3D electric and 
magnetic fields. Correspondingly, there is no limit on the initial 
current density distribution of the electron beam, nor its density 
distribution at any time during the simulation. 

Simulation results including beam current density, percent 
ripple and percent transmission will be presented, and the 
effects current, magnetic focusing strength and thermal 
velocities have on beam behavior will be demonstrated using 
3D movies showing the evolution of beam characteristics in 
time and space. Unlike typical beam optics models, this 3D 
model allows simulation of asymmetric designs such els non- 
circularly symmetric electrostatic or magnetic focusing as well 
as the inclusion of input/output couplers. 

1 The MAFIA Collaboration, MAFIA TS3 the 3 D-PIC Solver, December, 
1996. 
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♦ Current density 

♦ Tunnel length 
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♦ inclusion ot couplers 

♦ Input/Output sections 

♦ More accurate spent beam data for collector design 
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♦ Excessive growth in beam diameter leads to: 
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magnetic flux density B z on the PPM stack central axis 

versus z/L 



Longitudinal Magnetic Flux Density, Bz 



Percent ripple versus number of particles 



% Ripple 
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Thus, particle trajectories are emitted across beam cross section with random 
distribution of ejection angle 0 according to mentioned Gaussian probability 
function. 
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Percent ripple for various values of space charge 
parameter (3 at Brillouin focusing 



Percent Ripple -Comparison to 

Mendel et al. 




Percent ripple for 81mA beam with and without 
initial transverse velocities using various 

focusing strengths 




Percent transmission for 81mA beam for various values of 
flux density with and without initial transverse 

velocities 




Tunnel Length- Beam profiles/current density 
2*Brillouin with initial transverse velocities 
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♦ More accurate spent beam data for collector design 



